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SUMMARY 


Flat, single -Lay, skin stiffener panels with length-width ratios of O .96 
were tested at a Mach number of 3-0, at dynamic pressures ranging from 1,500 to 
5,000 pounds per square foot, and at stagnation temperatures from 300 F to 
660 ° F in order to investigate the effects of thermal stress and buckling on the 
flutter of such panels. The panel supporting structure allowed partial thermal 
expansion of the skins in both the longitudinal and lateral directions. Panel 
skin material and skin thickness were varied. 


A boundary faired through the experimental flutter points consisted of a 
flat-panel portion, a buckled-panel portion, and a transition point, at the 
intersection of the two boundaries, where a panel is most susceptible to flu er. 
The flutter region consisted of two fairly distinct sections, a large-amplitude 
flutter region and a small-amplitude flutter region. The results show that an 
increase in panel skin temperature (or increase in thermal stress) makes a flat 
(unbuckled) panel more susceptible to flutter. The flutter trend for buckled 
panels is reversed. Use of a modified temperature parameter, which approximately 
accounts for the effects of differential pressure and variations in panel skin 
material and skin thickness, reduced the scatter in the data which resulted when 
these effects were neglected. The results axe compared with an exact theory for 
clamped panels for the condition of zero midplane stress. In addition, a two- 
mode " tran stability" solution for clamped panels is compared with the experi- 
mentally determined transition point. 


INTRODUCTION 


Flutter characteristics of thermally stressed panels have become important 
with the advent of sustained supersonic flight and the aerodynamic heating 
associated with such flight conditions. However, very little experimental infor- 
mation exists concerning the flutter behavior of panels acted upon by compressive 
forces, particularly in the nonbuckled range. Some information on the effects of 
midplane compressive stress has been obtained experimentally through mechanical 
buckling (with some heating) of clamped plates having length-width ratios of 5 



or less. (See, for example, refs. 1 to 4 .) In addition, the effects of com 
pressive stress and buckling, induced by aerodynamic heating, on the flutter 
characteristics of multibay panels with length-width ratios of 10 have been 
reported in references 5 and 6. 


The experimental flutter trends presented in references 5 and 6 revealed 
that the flutter boundary consists of a flat-panel portion, a buckled-panel por- 
tion, and a transition point, at the intersection of the two boundaries, where a 
panel is most susceptible to flutter. The flat-panel boundary indicated that, 
for a given panel and aerodynamic conditions, an increase in panel temperature 
(or midplane compressive stress in the direction of air flow) made a flat panel 
more susceptible to flutter; the experimental flutter trends for flat panels were 
similar to trends predicted theoretically in references 7, 8, and 9. For the 
buckled-panel boundary the trends were reversed; once a panel is thermally 
buckled, any additional temperature rise will tend to stiffen the panel. 

The present investigation was conducted in the Langley 9 - by 6-foot thermal 
structures tunnel and was undertaken to determine the effects of compressive 
stress (induced by aerodynamic heating) and buckling on the flutter character- 
istics of panels with length-width ratios of O.96. Single-bay panels, 25.6 inches 
long and 26.6 inches wide, were tested at a Mach number of 5.0 at various dynamic 
pressures and stagnation temperatures. Panel skin material and skin thickness 
were varied. The panel- supporting structure was such that partial thermal expan- 
sion of the skin could occur in both the longitudinal and lateral directions and 
thus delay thermal buckling and permit more extensive study of the flat-panel 
flutter boundary . 


The flutter data obtained in this investigation are presented in tabular 
form. In addition, a modified temperature parameter, which approximately accounts 
for the effects of differential pressure and variations in panel skin material and 
skin thickness, is derived and used in conjunction with a flutter parameter to 
summarize the data as a flutter boundary which indicates the overall effects of 
compressive stress and buckling. The experimental data for flat, clamped panels 
are compared with results from an exact theory for the condition of zero midp lan e 
stress. In addition, a two-mode "transtability" solution for clamped panels is 
compared with the experimentally determined transition point . 
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spring constants per unit length (see fig. 12) 

panel length (longitudinal direction, parallel to air flow) 
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free-stream static pressure 
static pressure in bay behind panel 

differential pressure acting on panel skin. Pp - Poo 

dynamic pressure 
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initial temperature of panel 
stagnation temperature 
free-stream temperature 
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coordinate, parallel to airflow (see fig. 12) 
coordinate, perpendicular to air flow (see fig. 12) 
coefficient of thermal expansion of panel skin 
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specific weight of panel material 

average midplane stress in x-direction, positive in compression 
average midplane stress in y-direction, positive in compression 
thickness of panel skin 

modified temperature parameter, ±K S ^±a ^ - O.92 ^ 



TESTS 


Panels 

The panels^ which were 25-6 inches long and 26*6 inches wide; consisted of 
flat sheets of 0.072- or 0 . 08 l-inch-thick 2024--T3 aluminum alloy or 0 , 050 -inch- 
thick 17*7 PH stainless steel riveted to Z-section stiffeners along the longi- 
tudinal and lateral edges. The stiffeners were approximately 1.6 inches deep 
and were formed from 0 . 125 -inch- thick 2024-T3 aluminum alloy. Steel bars and 
angles were attached to the bottom of the stiffeners to provide support for 
mounting instrumentation. A rear-view photograph of a typical panel is shown in 
figure I, and pertinent panel construction details are given in figure 2. 


Test Apparatus 

Langley 9- by 6-foot thermal structures tunnel .- All tests were conducted in 
the Langley 9- by 6 -foot thermal structures tunnel, a Mach 5 intermittent blow- 
down facility exhausting to the atmosphere. A heat exchanger is preheated to 
provide stagnation temperatures up to 660° F and the stagnation pressure can be 
varied from 60 to 200 pounds per square inch absolute. (Additional details 
regarding the tunnel may be found in ref. 5») 

Panel holder and mounting arrangement .- The panels were mounted in a panel 
holder which extended vertically through the test section (fig. 5) • panel 

holder has a beveled half -wedge leading edge, flat sides, and a recess 29 inches 
wide, 50 inches high, and approximately 5*5 inches deep for accommodating test 
specimens. The recess is located on the nonbeveled side of the panel holder. 
Pneumatically operated sliding doors protect test specimens from aerodynamic 
buffeting and heating during tunnel starting and shutdown.. Aerodynamic fences 
prevent shock waves emanating from the doors from interfering with the airflow 
over the test specimen. (The flow conditions over the area of the recess are 
essentially free-stream conditions as determined from pressure surveys of a flat 
calibration panel (ref. 5)0 A vent-door arrangement on the side opposite the 
recess for the panel is used to control the pressure inside the chamber behind 
the test specimen. 

All panels were mounted flush with the flat surface of the panel holder. 

(See fig. 3.) The panels were attached to a mounting fixture (bolted to the 
panel holder) by means of 0 . 125 -inch-thick aluminum-alloy angles to provide 
support along the edges. (See fig. A.) This mounting arrangement, which allowed 
partial ther mal expansion of the panel skin in the longitudinal and lateral direc 
tions, was used in order to delay thermal buckling and thus permit more extensive 
study of the flat -panel boundary. 


Instrumentation 

Iron-constantan thermocouples, spotwelded to the panels at the l4 locations 
shown in figure 5, were used to measure panel temperatures. Inductance-type 
deflectometers were used to determine panel skin deflections. The deflectometers 

5 



were located approximately l/4 inch behind the panel skin at the four positions 
indicated in figure 5* In addition, high-speed l6-millimeter motion pictures 
provided supplementary data on panel behavior. Grid lines were painted on the 
panel skins for photographic purposes. 

Quick-response, strain-gage-type pressure transducers were used to measure 
static pressures in the tunnel and at various locations on the panel holder and 
in the chamber behind the panels . Tunnel stagnation pressures were obtained from 
static pressures measured in the settling chamber. Stagnation temperatures were 
measured by total temperature probes located in the test section. For each test 
all temperature and pressure data for both the test panels and the tunnel oper- 
ating conditions were recorded on magnetic tape. Deflection data were recorded 
on high-speed oscillographs . 


Test Procedure 

The tests were conducted at a Mach number of $. 0 , at dynamic pressures from 
1,500 pounds per square foot to 5,000 pounds per square foot, and at stagnation 
temperature s from 500° F to 66 0° F. The protective doors on the panel holder 
were opened after desired test conditions were established. The dynamic pressure 
was held constant during the first few seconds of all tests but was varied during 
the remainder of most tests. The usual procedure for varying the dynamic pres- 
sure was as follows: (a) if no flutter had occurred after a predetermined length 

of time, the dynamic pressure was increased in an attempt to initiate flutter; 

(b) if flutter had started and stopped, the dynamic pressure was increased in an 
attempt to restart flutter; (c) if the panel was fluttering near the end of a 
test, the dynamic pressure was decreased in an attempt to stop flutter. The 
differential pressure on the panels was controlled manually in order to keep the 
differential pressure as small as possible (usually less than 0.6 pound per 
square inch) . The stagnation temperature was essentially constant during most 
of the tests but decreased slightly during the latter portions of some tests. 

The protective doors were closed 5 seconds prior to tunnel shutdown. The dura- 
tion of test conditions varied between approximately 5 and 50 seconds. 


RESULTS AND DISCUSSION 


Seventeen tests were made in this investigation, fourteen on aluminum-alloy 
panels and three on steel panels . Flutter, initiated by an increase in panel 
skin temperature, dynamic pressure, or both, was observed in all tests. Pertinent 
data for all tests are given in table I. The data tabulated include the stagna- 
tion temperature T t , dynamic pressure q, differential pressure Zip, the skin- 

temperature increase ZiT, and the frequency f at the start of large -amplitude 
flutter . 
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Panel Temperatures 


At the beginning of a test the panel skin and supporting structure were 
essentially at the same temperature. Any temperature increase of the panels 
prior to opening the panel holder protective doors was usually insignificant. 
After a panel was exposed to the airstream, the skin temperature increased in a 
manner similar to the typical temperature histories shown in figure 6 (test l4) . 
The top curve represents the average of thermocouples 4 to 10 and l4j individual 
temperatures were within 10° F of the average value. The temperature histories 
for thermocouples 1, 2 , 3, 11, 12, and 13 indicate that there were appreciable 
lateral and longitudinal temperature gradients in the panel skin near the sup- 
porting structure and large temperature gradients in the supporting structure. 
However, these temperature variations were neglected in analyzing the test data 
and the average increase in temperature of thermocouples 4 to 10 and l4 was con- 
sidered to be the temperature increase Zff of the panel skin. For two tests 
(tests 9 and 13 ) all or most of the skin thermocouples were lost during testing. 
Hence, for these tests calculated temperatures, obtained by the same procedure 
as used in reference were used. Temperatures were calculated from the 
equation (ref. 10 ) 


-h(t-ti) 

T - - (^av - %)• ^ (D 

which neglects temperature variation through the skin, heat flow by conduction, 
and heat transfer by radiation. The aerodynamic heat-transfer coefficients were 
obtained from the turbulent -flow flat-plate theory presented in reference 11 by 
using initial free- stream flow conditions, a skin temperature equal to T^, and 
a reference length equal to the distance from the leading edge of the panel 
holder to the center of the panel. Temperature calculations, using adiabatic 
wall temperatures obtained in the usual manner ((T aw ) = q r (T t - tJ + T , 

where q r is the turbulent -flow recovery factor), gave skin temperatures higher 

than the measured temperatures. An arbitrary adjustment of the adiabatic wall 
temperature was therefore made "by using the relation 

MadJ = °-90(T aw ) calc (2) 


Use of this entirely empirical relationship in equation (l) resulted in fairly 
good agreement between measured and calculated temperatures. Figure 6 shows a 
comparison of measured and calculated temperatures for a typical test (test lk) . 


Flutter Parameters 

The flutter data obtained in this investigation are summarized in terms of 
a dimensionless flutter parameter and a dimensionless modified temperature 
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parameter. The flutter parameter 


PE/ 


1/5 


is proportional to the cube root of 


the primary panel flutter parameter given by theory (see, for example, ref. 7) 
and is the reciprocal of the parameter used in previous experimental investiga- 
tions. (See, for example, refs. 5 and 6 .) The dynamic pressure q, skin thick- 
ness t, and Young's modulus E were varied in this investigation. Because of 
the short duration of the tests, changes in material properties with temperature 
were assumed to be negligible . 


In this investigation the manual control of Ap resulted in differential 
pressure loadings of such magnitude that, when the data were plotted in terms of 

the parameter a AT^-^ 2 , as was done in reference 5 , there was considerable 


scatter in the data, and thus it was desirable to consider the effects of Ap. 
Therefore, the effects of midplane stress and buckling are indicated, in terms of 
the temperature rise AT and the differential pressure Ap, by the modified 
temperature parameter \|r where 


* = ±K S 



- 0.92 




( 5 ) 


In equation (3), which is derived in the appendix, the positive signs apply 
when a panel is flat inasmuch as \|r is then an indication of the stress ratio 

— ^ — and may be positive (compression) or negative (tension) depending on the 

x^D/l 2 

relative magnitudes of the temperature and differential pressure terms. However, 
when a panel is buckled, the negative signs apply inasmuch as the parameter \|r 
is then an indication of the buckle depth and both AT and Ap tend to increase 
the depth of buckle. (For this condition the panel is subjected to compressive 
stress and the negative sign in front of K s insures that \| r will be positive.) 

The effects of the lateral and longitudinal inplane restraint of the supporting 
structure on the midplane stress are approximated through the support factor K g . 

Values of K s are listed in table I and were obtained from equation ( A 6 ) ; the 
appropriate stiffness of the supporting structure was represented by spring con- 
stants kjj. and ky (per unit length) which were obtained from an approximate 

analysis of the bending stiffness of the panel supporting structure and support 
angles . 


Although theoretical investigations of the effects of midplane stress on 
panel flutter have used the ratio of midplane stress to Euler's buckling stress 

(for pinned ends) - — , preference for the parameter V for this investiga- 

1 r 2 D/l 2 

tion resulted from several considerations. The skin- stiffener type construction 
and mounting arrangement used for the panels of this investigation permit 
restrained inplane and rotational displacements at the edges of the panel due to 
both thermal expansion of the skin and differential pressure loading. The 
resulting vertical displacements of the panel skin would necessitate use of 
nonlinear large deflection theory in conjunction with rather complicated boundary 
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conditions in order to determine the midplane stresses correctly. Rather than 
resort to such calculations , the parameter \|/, which gives an indication of the 
N x 

stress ratio t-tt* was used to present the data, 

* 2 d !i 2 ’ 


Effects of Thermal Stress and Buckling on Panel Flutter 


Results from all tests are presented in figure 7 in terms of the flutter 
parameter and the modified temperature parameter. When a panel was flat, flutter 
occurred initially as a large -amplitude motion that was periodic but nonharmonic. 
The open symbols in figure 7 represent the start of this large -amplitude flutter 
resulting from increasing thermal stress under either constant or increasing 
dynamic pressure. In most tests, the large -amplitude motion changed to small- 
amplitude flutter or an intermittent small -amplitude oscillation as the panel 
skin temperature increased. The half -solid symbols represent this transition 
from large -amplitude flutter to small-amplitude flutter. When the small -amplitude 
flutter stopped as a result of further heating or a reduction in dynamic pressure, 
the panels appeared to be buckled. The solid symbols represent this cessation of 
small-amplitude flutter. When a panel was thermally buckled at the start of 
flutter, the ensuing motion was of the small- amplitude type. The open symbols 
with tick marks represent the start of small -amplitude flutter from the buckled 
condition. Thus, as can be seen from figure 7; the overall boundary, faired 
through the experimental flutter points, consists of a flat-panel portion, a 
buckled-panel portion, and a transition point at the intersection of the two 
boundaries. The flutter region, which is above the boundaries, consists of two 
fairly distinct sections, a large -amplitude flutter region and a small -amplitude 
flutter region. 


Flat-panel boundary and large -amplitude flutter region .- To the left of the 
transition point < 16) in figure 7* a panel is flat (^unbuckled) when flutter 
starts. For a given panel and aerodynamic conditions, an increase in compressive 
thermal stress will make a flat panel more susceptible to flutter as indicated by 
the downward slope of the flat -panel boundary. Thus, as the value of the tempera- 
ture parameter increases, the value of the flutter parameter required for flutter 
decreases until the transition point is reached. At the transition point, the 
value of the flutter parameter is a minimum ( 2 . 80 ) and is approximately 76 percent 
of the value for no stress (5*68). 


The flutter-start points (panel flat) for the steel panels fall within the 
scatter of the data for the aluminum-alloy panels, as shown in figure 7* Thus, 


it appears that the flutter parameter 


jj \ i/5 . 

PE j 


— and the modified temperature 

T 


parameter t|t adequately account for variations in panel skin material and thick- 
ness, at least for the flat (unbuckled) panel data for these tests. Since the 
Ap term in the parameter t|t greatly reduced the scatter that resulted when Ap 


was neglected (see appendix), use of the term 


f £)1 


2/3 


appears to be justi- 


fied as an approximation to the effects of Ap. Much of the remaining scatter in 
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the data probably results largely from the approximations used in calculating 
the effects of Ap and from the limitations imposed on i|f by neglecting certain 
factors such as edge rotation. 

In al l instances when a flat panel became -unstable due to an increase in 
thermal stress, the resulting flutter consisted of large -amplitude motion that 
was periodic but nonharmonic. The initial motion was relatively steady when 
flutter started and the fundamental flutter frequency varied from 70 cycles per 
second to l40 cycles per second. The basic mode appeared to consist of a single 
half-wave in both the lateral and longitudinal directions. However, some higher 
frequency, lower amplitude motion was nearly always present as indicated in fig- 
ure 8(a). (it should be noted that all deflectometers were located off center of 
the panel because the amplitude of the motion was so large that the panel skin 
could come in contact with deflectometers located near the center of the panel.) 

As a test progressed, the motion became more irregular as the higher frequency 
portion of the motion increased (fig. 8(b)). The irregular motion either con- 
tinued to the end of the test or until there occurred a fairly distinct decrease 
in flutter amplitude (fig. 8(c)), which was usually accompanied by a change to a 
flutter mode consisting of several half-waves in one or both directions. The 
transition to small-amplitude flutter is indicated by the half-solid symbols in 
figure 7. As can be seen from the figure, the points representing the transition 
to small-amplitude flutter for the steel panels are in poor agreement with the 
data for the aluminum- alloy panels. This condition suggests that the parameter \|/ 
does not account for variations in panel skin material for buckled panels . 

Buckled-panel boundary and small-amplitude flutter region .- The solid 
boundary to the right of the transition point U > l6, figT 7 ) represents the 
buckled-panel flutter boundary, below which the panel is buckled and stable. As 
can be seen from the figure, the points representing the start of small-amplitude 
flutter (open symbols with tick marks) agree fairly well with the boundary estab- 
lished by the small- amplitude flutter stop points (solid symbols). The slope for 
the buckled boundary is reversed from that of the flat-panel boundary. This same 
reversal in trend has been shown experimentally for panels with length-width 
ratios of 10 (refs. 5 and 6) and theoretically by Frail ch (contained in ref. 12) 
for a simply supported panel with length-width ratio of 1. The reversal in trend 
is attributed to the panel post-buckling behavior. Once a panel becomes thermally 
buckled, any additional temperature rise will increase the depth of buckle which 
tends to stiffen the panel (due to lateral curvature) as shown both theoretically 
and experimentally in terms of panel natural frequencies in references 13 and l4. 
In the stable buckled-panel region, and well to the right of the transition point, 
the buckling mode had either 1 or 3 half-waves in the longitudinal direction and 
3 half-waves in the lateral direction (fig. 9) • 

The significance of the dashed boundary between the large -amplitude and 
small-amplitude flutter regions and the variations of panel flutter behavior 
within the s mall -amplitude flutter region are not clearly understood . To the 
left of the dashed boundary, the flutter mode always appeared to consist of a 
single half -wave in both the longitudinal and lateral directions, which corre- 
sponds to the fundamental buckling mode for square panels (for Oyja x = 1.0). To 
the right of the dashed boundary the flutter amplitude was small and the motion 
inconsistent. For several tests, flutter was observed to stop and then restart 
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(occasionally several times for a given test) within the small -amplitude flutter 
region* When flutter stopped, the panel appeared to be buckled. However, both 
the flutter and buckling modes were observed, at different times, to have 1, 2, 
or 3 half-waves in both the longitudinal and lateral directions. Usually the 
small-amplitude flutter consisted of low- frequency motion near the leading edge 
and high-frequency motion near the trailing edge but the motion did not appear to 
be entirely superimposed on the buckles. However, in some instances the motion 
consisted of very low- amplitude, high-frequency oscillations superimposed on the 
buckles or on a low-frequency oscillation of the buckles. In other tests the 
motion degenerated to a low-amplitude, intermittent oscillation that was aperi- 
odic. In addition, short bursts of large -amplitude flutter, usually accompanied 
by an abrupt change in mode shape, were observed in the small-amplitude flutter 
region. In many instances, several of the different types of motion described 
above were observed in the small -amplitude flutter region during a given test. 

These variations in flutter behavior within the small -amplitude flutter 
region appeared to be inconsistent and no definite trends could be established 
from the data available. Therefore, these variations in flutter behavior were not 
indicated in figure 7> and only those points used to establish the boundaries 
that define the small- amplitude flutter region are shown. One factor which prob- 
ably contributed to the variations in panel-flutter behavior was the variations in 
buckling mode; variations in buckled mode shape have been shown experimentally 
(ref. l) to affect the flutter characteristics of panels. Stein (refs. 15 and l6) 
has shown that increases in loading of buckled elastic structures (including 
rectangular plates) can result in changes in the buckled mode shape. Thus, the 
variations in panel flutter behavior suggest the possibility of additional 
buckled-panel boundaries within the small -amplitude flutter region. Moreover, 
the additional data and more thorough knowledge of the buckled mode shapes 
required to define these boundaries might also alter the boundaries shown in fig- 
ure 7, which were arbitrarily faired through the data that defined the bounds of 
the small -amplitude flutter region. 

That the boundary between the large -amplitude and small-amplitude flutter 
regions can be crossed from either side is suggested by the results of one test, 
as shown in figure 10. This figure shows the flutter boundaries of figure 7 an d 
the variation of the flutter and modified temperature parameters during test l4 
wherein the panel became thermally buckled prior to flutter. The panel, which was 
initially flat and stable, buckled in the vicinity of the transition point. (The 
buckle depth was so small that the number of half-waves could not be definitely 
ascertained.) The dynamic pressure was then increased and small- amplitude flutter 
started as indicated by the open symbol with tick mark. As the dynamic pressure 
continued to increase there occurred (near the dashed boundary) a transition from 
small-amplitude flutter to large-amplitude flutter and back to small- amplitude 
flutter, as indicated in figure 10 . (The burst of large -amplitude flutter was not 
shown in fig. 7j since all other flutter points establishing the dashed boundary 
represented transition from large -amplitude flutter to small -amplitude flutter.) 

One panel was damaged while undergoing small-amplitude flutter (test 13) . 

The damage, which occurred approximately 6 seconds after flutter began, started 
as a tear along the inner rivet line at the trailing edge. This was quickly fol- 
lowed by several tears parallel to the airflow, which resulted in shredding of 
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the panel skin (fig. 11 ) . The panel had previously been subjected to approxi- 
mately 17 seconds of flutter (test 9) • 


Comparison with theory for conditions of zero stress and onset of buckling . - 
A theoretical flutter boundary, based on an exact solution for flat, finite, 
clamped panels subjected to midplane stress, is presented in reference 8 in terms 

2 qZ 5 — 

of the parameters — and A* However, the boundary of reference 8 extends 


only to A = -1, whereas for__the condition of zero stress for the panels of this 
investigation the value of A is -2.46. Thus it was necessary to make additional 
numerical calculations in order to extend the boundary of reference 8. The theo- 
retical value of the flutter parameter so obtained for zero midplane stress was 
5.40 ( see fig. 7) which is 92 percent of the experimental value of 3.68 (in 
obtaining the theoretical value, the Mac h number M, which appears in the param- 
eter of ref. 8, was replaced by \/m 2 - 1 which is applicable for a Mach number 
of 3.0) . 


Isaacs (ref. 17) introduced the concept of the "transtability" speed, a 
speed, calculated from purely static considerations, that constitutes the flutter 
speed of a panel at the onset of buckling. The transtability concept has been 
used in two mode analyses to determine the effects of the onset of buckling (which 
defines the transition point) on simply supported panels (ref. 7) and clamped 
panels (ref. l8) . Reference 7 indicates that the value of the flutter parameter 
at the transition 'point depends on the stress ratio dyycr x ; this ratio was esti- 
mated to be 1.0 for the panels of this investigation. The transition point for a 
square, clamped panel (<jy/cx x = 1.0) has been determined by the "transtability” 

method in reference l8 and found to be 2.20 (see fig. *]) > which is approximately 
79 percent of the experimental value of 2.80. The fact that the agreement between 
theory and experiment is not as good for the onset of buckling as for the condi- 
tion of zero stress may be due to several factors. As the stress a x increases 

the effect of the differences in the actual and theoretical panel boundary condi- 
tions becomes more pronounced since the edge rotation induces vertical displace- 
ment (or curvature) of the panel skin. In addition, a two-mode analysis would be 
on the conservative side of the exact solution (refs. 7 and 8); in this case, a 
two-mode analysis increases the discrepancy between theory and experiment. In 
any event, it appears that the theory is conservative and is In fair agreement 
with experiment, at least for the panels of this investigation. 


CONCLUSIONS 


Flat, single-bay, skin- stiffener panels with length- width ratios of O.96 
were tested in the Langley 9- by 6- foot thermal structures tunnel to investigate 
the effects of thermal stress and buckling on the flutter characteristics of such 
panels. The tests were conducted at a Mach number of 3*0, at dynamic pressures 
ranging from 1,500 pounds per square foot to 5>000 pounds per square foot, and at 
stagnation temperatures from 300° F to 660° F. The panel supporting structure so 
restrained the panels that partial thermal expansion of the skins could occur in 
both the longitudinal and lateral directions. The ratio of the lateral to 
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longitudinal midplane stress in the panel skins was estimated to be approximately 
1.0. The panel skins were 0.072- or 0 .08l-inch- thick aluminum alloy or 0.050- 
inch-thick stainless steel. A dimensionless modified temperature parameter, 
which approximately accounts for the effects of differential pressure and_ varia- 
tions in panel skin material and skin thickness, was derived and used to indica e 
the flutter characteristics of thermally stressed panels. The tests revealed the 
following: 

1. The overall flutter boundary consisted of a flat-panel portion, a buckled- 
panel portion, and a transition point between the two. The flutter trends indi- 
cated by the overall boundary were similar to the experimental trends obtained 
previously for panels with length-width ratios of 10 and to theoretical results 
obtained for simply supported panels with a length-width ratio of 1.0. 

2. At the transition point, the value of the flutter parameter was a minimum 
and was approximately 76 percent of the value required for flutter of an unheated 

panel. 

5 . The flutter region consisted of two relatively distinct sections, a large- 
amplitude flutter region and a small-amplitude flutter region. Panel behavior m 
the large-amplitude flutter region consisted of periodic but nonharmonic motion. 
Panel flutter behavior in the small -amplitude flutter region was inconsistent. 

The variations in panel flutter behavior (in the small-amplitude flutter region) 
suggest the possibility of additional buckled-panel flutter boundaries within, 
this region. For small- amplitude flutter, the panels usually appeared to be in a 
huckled condition. 

4. The flat -panel (unbuckled) flutter data for the steel panels fell within 
the scatter of the data for the aluminum- alloy panels and indicated. that the 
parameters used to summarize the data adequately account for variations in panel 
skin material and thickness. 

5. Use of a differential pressure term in the modified temperature parameter 
greatly reduced the scatter in the data that resulted when the differential pres- 
sure Ap was neglected and thus the use of this term appears to be justified as 
an approximation to the effects of Ap. 

6. Exact theoretical results for finite, clamped panels were in fair agree- 
ment with the experimental results for the condition of zero midplane stress, 
two-mode "transtability" solution for clamped panels (for the onset of buckling) 
was in fair agreement with experiment at the transition point. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 25, 19^2. 
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APPENDIX 


DERIVATION OF MODIFIED TEMPERATURE PARAMETER 


As noted in the section entitled "Flutter Parameters/' a dimensionless 
temperature parameter (modified to include the influence of differential pressure) 
is used to show the effects of both aerodynamic heating and differential pressure 
on the flutter data. The analysis herein uses stress-displacement relationships 
to relate panel midplane stress to temperature (see, for example, ref. 8) but 

employs a semiempirical approach for the inclusion of the effects of differential 
pressure. 


Consider the idealized rectangular panel shown in figure 12. The panel 
edges are not permitted to deflect out of their original plane and longitudinal 
and lateral expansion is restrained by springs (with spring constants k x and k y 

per unit length, respectively) which represent the appropriate stiffness of the 
supporting structure. The panel skin is assumed to be subjected to a uniform 
temperature increase AT, thermal expansion of the supporting structure is 
assumed to be negligible, and a uniform differential pressure Ap acts on the 
panel so as to produce vertical displacement 5. Under these conditions, prior to 
buckling, compatibility in the x- and y-directions requires that 


P X T 
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(a AT + 




n 
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(Ala) 
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— fa AT + tier 
2y E x 
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yields 


Solving equations (Al) for the average midplane stresses a x and a 
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where 


and 


e x 


i 



0y 
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+ 2E t 

Vky 


(A3a) 


(A5b) 


If it is assumed that the panel is square (l = w) and that k x = ky, then 
p x = Py, and 7 X = 7y and equation (A2a) may be written 


= 




(A4) 


Substitution of equation (A4) into the frequently used dimensionless stress ratio 
(the ratio of midplane stress in the direction of airflow to Euler's buckling 
stress for pinned ends) yields 


where 


N. 


x 


rt 2 D /l 2 


= K c 



K s = 


12(1 - p 2 ) 

« 2 (Px " d) 


(A5) 


(A6) 


An approximation to the effects of Ap can be obtained by considering only 
the direct effect of membrane stresses induced by Ap; this approach neglects 
such effects as the interaction of induced thermal stresses and vertical displace- 
ments. Thus, if it is assumed that the deflected shape is a function of Ap 
only, then, since / x is a function of the deflected shape, equation (A5) may be 
written 


N x 

it 2 D h 2 


= K„ 


a AT 


- C 1 g(Ap) 


(A7) 


In this expression g(Ap) is a nondimensional function of Ap, upon which 

depends when AT = 0, and C-i is a proportionality constant that can be 

*2D/l2 
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determined from the experimental data. The form of the function g(Ap) can he 
obtained from an analysis of the stress in a square membrance subjected to a 
uniformly distributed normal loading. (See ref. 19.) When this stress is sub- 
stituted in the stress ratio, the stress ratio is seen to be proportional to the 


dimensionless parameter 



and thus equation (A7) may be expressed as 


rt^/z 2 




(A8) 


The constant Cp, obtained from the flat-panel flutter data, was found to be 
0.92. 


The results of the present invest! gat ion are presented in terms of the two 
parameters on the right-hand side of equation (A8) rather than in terms of the 
stress ratio itself in order to emphasize the approximations used in evaluating 
the effects of both aerodynamic heating and differential pressure and in order 
not to imply that the midplane stresses were accurately evaluated. Thus, equa- 
tion (A8) becomes 



The positive signs apply when a panel is flat as t is then an indication of the 
N x , 

stress ratio and may be positive (compression) or negative (tension) 

jt 2 d/ Z 2 

depending on the relative magnitudes of the temperature and differential pressure 
terms. When a panel is buckled, the negative signs apply as the parameter \|r is 
then an indication of the buckle depth and both AT and Ap tend to increase 
the depth of buckle . (For this condition the panel is subjected to compressive 
stress and the negative sign in front of K s insures that \|/ will be positive.) 

The parameter i|f was arbitrarily chosen to present buckled-panel data in order 
to avoid an additional and more complicated analysis; moreover, use of the param- 
eter ijr tended to reduce the scatter in the data for the buckled-panel 
boundaries . 


The constant Cj_ was determined as follows. All flutter start points 

) l/ 3 

— 

T 

and K s a AT f Aj 2 . An arbitrary boundary was faired through the points for which 
Ap_ was small (Ap 5 0.08) as shown in figure 13. Then the values of 

(^) J , representing the differences between the flutter points (to the 


K s a AT( 
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right of the boundary) and the arbitrary boundary, were plotted against the 

fAp/A^l 2 / ^ 

appropriate values of the differential pressure parameter 1 — * — * 1 


[fm 


Since 


K^a AT( — ] I is an indication of the contribution of Ap to the stress ratio 


i — 

n x r 

A L Ksa 


AT^- 


might be expected to vary approximately linearly with 


&p/i 
E \T 



2/5 


. The method of least squares was used to obtain a straight line 


through the data; the slope of this line ( 0 . 90 ) gave a first approximation to the 
constant C^. The experimental flutter points were replotted by using the param- 


eters 


1/5 

pE/ T 


!V " 1 and K c 


a AT^i 


boundary. Then values of A 

i2 


K s a AT[1- 


- 0.90] 

n 2 


Ap (j}^ 


2/5 


to obtain a new 


the values of K, 


,a 


, representing the differences between 
at the flutter points and the values of 


Ks - o. 9 ojf(i) 

2/5 


2/51 


at the new boundary, were plotted against 
to obtain a final value for of O .92 as shown in figure l4. 


Use of the parameter \|r (with the aforementioned change in sign) tended to 
reduce the scatter in the buckled-panel data to the extent that reasonable 
flutter boundaries defining the small-amplitude flutter region could be estab- 
lished; thus, this procedure was not performed for the buckled-panel data and the 
value of of O .92 was used in the parameter \|r to present all flutter data. 
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TABLE I.- PANEL FLUTTER DATA 
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pressure increasing 
pressure decreasing 






ro 


L-6l-5T 1 +5.i 

Figure 1.- Back view of typical panel (with support angles attached) 





Panel symmetrical about center line 



Section A— A Section B— B 

Figure 2.- Panel construction details. (All dimensions are in inches.) 
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L-61-3962.1 

Figure 3.- Panel mounted in panel holder in test section as viewed from 
upstream. Panel holder protective doors are in open position. 
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Figure 4 .- Panel mounting detail. Typical for all edges. 








Measured 

Calculated 



Protective doors open 

<£>Average of thermocouples 4 through 10 and 14 
□ Average of thermocouples 3 and II 
A Average of thermocouples 2 and 12 
OAverage of thermocouples I and 13 


Figure 6.- Measured and calculated panel temperatures for test l4 






(t>) Irregular large -amplitude flutter. 



(c) Transition to small -amplitude flutter under 
constant dynamic pressure. 

Figure 8.- Sample deflectometer records. 
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L- 62- 2122 

Figure 9 .- Panel buckled mode shape at cessation of small-amplitude flutter (test 8 ). 


Large— amp I i tude flutter start 



Figure 10.- Variation of flutter and modified temperature parameter during test in which 
the panel was thermally buckled prior to flutter. Positive signs in the expression 
for \|r apply when panel is flat, negative signs when panel is buckled. (See appendix. 




L-61-57^-1 

Figure 11.- Damage to panel (downstream half) with 0.072-inch-thick 
aluminum-alloy skin after 6 seconds of flutter (test 13). Panel 
had previously fluttered for 17 seconds (test 9). 
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Figure 12.- Idealized panel configuration -where panel is subjected to 

uniform temperature increase AT and uniform differential pressure Ap. 
Spring constants k^ and ky represent effective stiffness of unheated 
supporting structure . 
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Figure 13.- Approximate flat -panel flutter boundary where datp 
unmodified for effects of differential pressure. 
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Figure 1^.- Linear variation of A 
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